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In this work we performed a detailed investigation about ordering in Ba3CaNb20g 
perovskite. The sintering temperature and time were changed to obtain sam- 
ples with different ordering. The order parameters were probed by Raman 
spectroscopy based on a partial disordered model. To use the partial disor- 
dered model correctly we performed ab initio calculations in BaaCaM^Og to 
assign the optical phonons. The results showed that sintering temperature 

■"sj- ! improves order while sintering time is not so efficient to promote order. 
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1. Introduction 

Complex triple perovskites with the general formula AsB'B'^Og, where 
A is an alkaline earth metal, B' and B" are ions with valencies 2+ and 5+, 
have recently attracted much interest as high-temperature proton conductors 
[1, 2, 3, 4, 5, 6, 7] and as dielectric resonators [8, 9, 10, 11, 12, 13, 14]. 
Particularly, the nonstoichiometric series Ba3Cai +x Nb2xO(9_3 :r )/2 has been 
extensively investigated due to its high proton conductivity and chemical 
stability [1, 4, 5, 6, 7]. In such oxides, the formation of oxygen vacancies is 
achieved by increasing the B' / B" ratio. The vacancies can be filled with OH 
species upon annealing samples in water vapour at high temperatures and 
even OH - concentration is small, the ionic conductivity can be dominated 
by proton transport [15], enabling it for application in fuel cells. 

Despite oxygen vacancies can be highly mobile in simple ABO3 perovskites, [16] 
in complex perovskites as BaaCaNh^Og the vacancy mobility tends to be 
strongly influenced by the B 1 / B" ratio, valencies and ordering. Substitu- 
tions in A and B sites of ABX 3 perovskites are responsible for improving 
a large part of the applications of these materials, bringing new physical 
and chemical properties. [8] However, the substitution by more than one ion, 
mainly in B site, can lead to structures which are disordered or ordered, ac- 
cording the new ions distribution into the structure. This substitution leads 
to new crystalline structures when there is a ordered distribution. B-site 
cation-ordered triple perovskites whose chemical formula is A3B'B 2 'Og are a 
particular family obtained by B-site substitution . 

From the structural viewpoint, in the triple B-site ordered substitution 
(sometimes called 1:2 order), two main structures can be formed. In the First, 
an ordered B'/B" substitution into B site can lead to a trigonal structure 
which belongs to the space group P3ml. In this structure the cations B' and 
B" are alternately distributed in {111} planes in the form -B'B" B" B 1 B" B" B '■ 
■-. Another possible structure for this composition is the so-called hexagonal 
6H structure that belongs to the PQ^/mmc space group, in which the B'06 
and B"06 octahedra share edges. [7, 12] 

However, partial disordered B'/B" distribution in both structure can oc- 
cur, with part of the B' cations occupying the B" cation sites and vice versa. 
Partial disorder in perovskites plays an essential role because they strongly 
influence their physical properties. For example, order influences phonons 
and dielectric constant, implying in consequences for applications of 1:2 per- 
ovskites in wireless communication system. [17, 18, 19, 20, 21, 22, 11, 9, 23] 



Disorder is particularly interesting in Ba 3 CaNb 2 09 perovskites because it 
favors a higher conductivity in these protonic conductors based on nonstoi- 
chiometry Ba compounds, which always have higher conductivity and lower 
activation energy than the corresponding Sr compounds. [2] Due to the im- 
portance of ordering in BasCaM^Og (BCN) ceramics, mainly because its 
conductivity properties, in this work we proposed a systematically investiga- 
tion of the order achievement in BCN complex perovskite. 

2. Experimental Procedures 

BCN samples were synthesized by polymeric precursor method using bar- 
ium nitrate (Ba(N03)2, Sigma Aldrich) , calcium citrate tetrahydrate (Ca3(C6H 5 07)2-4H 2 0, 
Ecibra) and ammonium complex of niobium (NH 4 (NbO(C204)2(H 2 0)2) -3H 2 0, 
CBMM) as sources of metals [24]. Barium polymeric precursor was obtained 
by mixing aqueous solutions of barium nitrate and citric acid (CeHgO/- H2O, 
Proquimico) in a molar ratio of 1:3 metal-citric acid, keeping stirring between 
60 and 70 °C . Finally, ethylene glycol (HOCH 2 CH 2 OH, Merck) was added 
to metal/Citric aqueous solution in a mass ratio of 1:1 in relation to citric 
acid. Calcium polymeric precursor was obtained following the same proce- 
dure for barium. To obtain niobium precursor we firstly precipitated niobium 
oxi-hydroxide stirring an aqueous solution of ammonium complex of niobium 
until pH of 9 in a thermal bath at °C . Niobium hydroxide (Nb(OH) 5 ) was 
separated from oxalate ions using distilled water at 40-50 °C under vacuum 
filtering. Finally, citric acid and ethylene glycol were added to the Niobium 
hydroxide aqueous solution. The pH of the polymeric precursors was kept 
at the same values to avoid precipitations when mixing them. To determine 
the necessary precursor weight to achieve the correct metal stoichiometry to 
form BCN perovskite we used gravimetric analysis at 900 °C for 1 h. We 
mixed the three precursors and heated the mixture between 80 and 90 °C to 
form the polyester resin, which had high viscosity and glassy. The resin was 
annealed at 400 °C for 2h. This heat treatment converted the resin in a black 
porous powder. This powder was lightly grounded using an agate mortar. 
After it was calcined at 1300 °C for 2h to obtain BCN sample. This sample 
is our reference to investigate the ordering process and it will be referenced 
as start sample along this work. From this reference sample, several pellets 
were sintered at different temperatures and time of sintering to investigate 
the ordering phenomena. 

The crystalline structure of the samples were probed by powder X-ray 



diffraction (XRD - Bruker D8 Advance). We performed a continuous scan- 
ning mode using Cu-K al radiation (40 kV, 40 mA) over a 29 range between 
10° and 100° (0.02°/step with 8 seconds/step). The powder XRD pattern 
was compared with data from ICSD (Inorganic Crystal Structure Database, 
FIZ Karlsruhe and NIST) International diffraction database (ICSD#162758). 
[25] The structure was refined using the DBWS9807 free software. [26, 27] 

The Raman spectra of the samples were acquired at room temperature in 
an iHR550 Horiba scientific spectrometer coupled to an Olympus microscope 
model BX-41. A He-Ne laser (632.8 nm, 10 mW) was use to excite the spectra 
that were collected in an air-cooled Synapse CCD detector. The spectral 
resolution was kept lower than 2 cm" 1 using an 1800 grooves/mm grating in 
the spectrometer. All spectra were acquired in a backscattered configuration. 

3. Computational method details 

The first-principle calculations were performed through the Cambridge 
Serial Total Energy Package (CASTEP) software package [28]. Density func- 
tional theory (DFT) [29, 30] was chosen to model the material consider- 
ing generalized gradient approximation (GGA) as exchange-correlation func- 
tional. This approximation was improved by the Perdew-Burke-Ernzerhof 
(PBE) parameterization [31]. The PBE functional leads to results close to 
those obtained using PW91 functional [32] within TS method [33] for the 
dispersion correction scheme (DFT+D) to describe van der Waals interac- 
tions. We also adopted pseudopotentials to replace the core electrons in each 
atomic species. Specifically, norm-conserved pseudopotentials [34] were used. 
These pseudopotentials were generated using the OPIUM code [35], follow- 
ing the same scheme of previous works [36, 37, 38]. The electronic valence 
configurations for each atomic species were: Ba— 5s 2 5p 6 6s 2 , Ca— 3s 2 3p 6 4s 2 , 
Nb-4d 4 5s\ and 0-2s 2 2p 4 . A Monkhorst-Pack [39] 3x3x2 sampling was 
used to evaluate all integrals in the reciprocal space. This sampling is enough 
to give a well converged electronic structure. 

We employed BFGS minimizer [40] to optimize the unit cell. In this 
minimization algorithm, a starting Hessian is recursively updated. For each 
self-consistent field step, the electronic minimization parameters were: Total 
energy/atom convergence tolerance of 0.5 x 10~ 6 eV, a maximum energy 
eigenvalues threshold of 0.1442 x 10~ 6 eV, and a convergence window of 3 
cycles. A plane-wave basis set was adopted to represent the Kohn-Sham 
orbitals, with cutoff energy chosen of 880 eV. This value was obtained after 



convergence studies. The quality of this basis set was kept fixed as the unit 
cell volume varies during geometry optimization. 

The vibrational properties of BCN were calculated performing density 
functional perturbation theory (DFPT) calculations or linear response for- 
malism [41]. The linear response provides an analytical way of computing the 
second derivative of the total energy with respect to a given perturbation. 
Depending on the nature of this perturbation, a number of properties can be 
calculated. A perturbation in ionic positions gives the dynamical matrix and 
phonons; in the presence of magnetic field it gives the NMR response; in unit 
cell vectors changes it gives the elastic constants; in an electric field presence 
it gives the dielectric response, etc. The infrared absorption intensities are 
described in terms of a dynamical matrix and Born effective charges (also 
known as atomic polarizability tensors, ATP) [41] and can be obtained by 
calculating the phonons at the T point (k — 0). The structure used for vibra- 
tional calculations was the fully ordered trigonal BCN structure optimized 
through the GGA-PBE approximation. The geometry optimization criteria 
were more rigorous than the criteria used for LDA-DFT calculations, for 
example. The convergence thresholds used were: Total energy convergence 
tolerance smaller than 5 x 10 -6 eV/atom, maximum ionic force smaller than 
10~ 2 eV/A, maximum ionic displacement tolerance of 5 x 10 -4 A, and max- 
imum stress component smaller than 2 x 10~ 2 GPa. For the self-consistent 
field calculations, the convergence criteria considered a total energy per atom 
variation smaller than 5 x 10~ 7 eV, and electronic energy eigenvalue variation 
smaller than 0.1442 x 1(T 6 eV. 

4. Results and discussion 

Figure 1 shows the powder X-ray diffraction pattern obtained for the 
start sample, which was sintered at 1300°C for 2h. The reflection planes 
were indexed according to the trigonal structure which belongs to the space 
group P3ml. The refinement parameters were summarized in Table 1. In 
this structure the cations Ca 2+ and Nb 5+ are alternately distributed in {111} 
planes in the form • • -Ca — Nb — Nb — Ca — Nb — Nb ■ ■■. This distribution 
produces characteristic reflections of the trigonal superstructure. [42] When 
Ca and Nb are randomly distributed, and consequently the structure is disor- 
dered, the powder pattern does not exhibit such reflections, inducing changes 
in the plane indexing. In this case, a cubic structure that belongs to the space 
group Pm3m is usually observed. 
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Figure 1: Powder X-ray diffraction pattern of BCN ceramic sintered for 1300 °C /2h. The 
solid line is the fitting using the Rictvcld method and difference between the experimental 
and calculated patterns. The indices of planes are labeled in picture. 



Space group; Z 

Lattice parameters, A 

Temperature, K 

Density (calculated), g/cm 3 

Volume, A 3 

A, A 

Profile function 

Cagglioti parameters (U, V, W) 0.2056(2), -0.0496(8), 0.0549(6) 

Background function 5th order polynomial in 26 

Rp, Rwp, Rexp, % 7.63, 10.05, 5.90 



P3ml (N° 164); 3 

a=5.8921(8), c=7.2358(8) 

298(1) 

5.965 

217.56 

1.54056 

Pseudo-Voigt 



Table 1: Data collection and Refinement Details for BCN sample sintered for 1300°C for 
2h. 
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Figure 2: Raman spectra of the BCN ceramics sintered for several values of sintering 
temperature and sintering time. 



To investigate the formation of ordered BCN structure we performed 
Raman scattering measurements in BCN samples sintered for several values 
of sintering temperature and sintering time. The observed Raman spectra are 
shown in Figure 2. To classify and identify the phonons observed we used a 
group theoretical analysis. Thus, we determined the symmetry and number 
of phonons in the full ordered trigonal structure based on the irreducible 
representation of the group factor m3m. [43] In this structure the primitive 
unit cell has 15 atoms, where one barium ion is at the la site and two others 
are at the 2d site, the calcium ion is at 16 site and the two niobium ions are 
at 2d site, six oxygen ions are at 6i site and three are at the 3e site. [25, 44] 
Therefore, based on this ion occupation, the expected phonons in BCN are 
shown in Table 2. However, nine Raman-active modes (4t4i 3 © 5E g ) are 
expected in BCN at room temperature. 

In another way, occupational disorder can occur in two different ways: It 
implies in extra optical-active modes due to the different sites occupied by 
the ions (two-phonon behavior), [45] or it implies in local symmetry lower- 
ing (one phonon behavior) , [46] changing the Raman and infrared spectrum. 
In this work we assumed the first case, that is the most common in 1:2 
perovskites.[47, 45] Thus, to index the BCN crystalline structure according 
to the trigonal structure implies in an order-disorder model in which the trig- 
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Table 2: Distribution of modes in the perovskite crystalline structure belonging to the 
trigonal space group P3ml. 



Ion Site Symmetry Distribution of modes 

Ca2 2d 3m A lg © A 2u ®E g ®E u 
Nb2 lb -3m A 2u @E u 

Table 3: Distribution of modes for additional sites in the order-disorder model for the 1:2 
pcrovskitc compounds partially ordered. 



onal perovskite structure is not fully ordered, although it still has a trigonal 
unit cell. [48, 47, 14, 49] In this case, there are part of calcium ion occupy- 
ing the site 2d and part of niobium ion at the site lb. [50, 51] Thus, in this 
order-disorder model, the group theoretical analysis can be used to predict 
the extra optical-active phonons based on extra site occupation by Ca and 
Nb ions, as showed in Table 3. Therefore, considering the disorder, there are 
now eleven Raman-active phonons in the partially ordered BCN. 

From the Figure 2 (See bottom spectra in Figures 2(a) and 2(b)) we ob- 
served fourteen vibrational modes in the start sample, confirming a partially 
ordered 1:2 perovskite structure [47, 45] (see Figure 2b). Clearly, BCN Ra- 
man spectra change significatively when the sintering temperature and time 
are increased with both playing an essential role in ordering process. 

The most ordered sample is achieved when sintered at 1600 °C for 2h, 
whose spectral deconvolution is shown in Figure 4 and summarized in Table 
4. 

To identify the Raman-active phonons in BCN and check the partial disor- 
dered model, we performed ab initio calculations of the vibrational properties 
of the full ordered BCN trigonal structure. The good reliability of the method 
employed can be seen comparing the experimental[25] and calculated struc- 
tural data given in Table 5. The calculated phonons are summarized in Table 
4 together with the experimental phonons. The infrared-active phonons are 
also shown for completeness. 

The occurrence of the peaks n° 10 and 12 is correlated to the partial 
ordering adopted by the BCN confirming the disorder model adopted. This 
behavior can be explained by the Nb-0 bond length decreases as consequence 
the Nb 5+ substitution by the Ca 2+ , resulting in a displacement of the Ai s 
(821 cm -1 ) mode to lower wavenumber. [44] Also, a simple harmonic model 

to these vibrations, based on the charge/mass relation ? Ca ,™ u/ 2 give us the 
value 0.96, showing that the A\ g mode due to the extra-site Ca 2+ occupation 
should be observed 0.96 times that A\ g band associate to the Nb 5+ correct 
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Figure 3: Deconvolution of the Raman spectra of the BCN pellet sintered at 1600°C /2h. 
Collected data (o) and calculated Lorentzian (-) curve by fitting. 
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Tabic 4: Raman- and IR-active modes of trigonal BCN structure with respective attribu- 
tion. A comparison between experimental (EXP.) and calculated (CAL.) data is shown to 
Raman modes. a FLB: Floating Base Line; 6 DAM: Defect activated IR-activc modes. 
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occupation. [47] The calculated ratios, between the peaks n° 12 and 13 (0.92) 
and n° 10 and 11 (0.93) (see Table 4), show this assumption is correct. 

The evaluation of the ordering kinetics of the BCN structure was followed 
according the evolution of the A lg modes under the sintering temperature 
and time changes. A detailed spectral deconvolution of these modes are 
shown in Figure 4. 

Clearly, the intensities of the peak n° 10 and 12 are proportional to the 
percentage of the extra-site Ca 2+ and Nb 5+ ions. Using the modes n° 12 
and 13 we can estimate the ordering degree setting a ratio designated ^ca,Nb 
expressed by 

*Ca,N b = j^r- (1) 

J 12 + J 13 

When \I/jv6 is equal to one, all Nb ions are in correct site. In this case ^ca is 
null. Thus, the fully ordered BCN structure occurs when ^ca = h 2 = and 
^Nb = Irs = 1- Figure 5 shows the evolution of the ratio ^ca,Nb under the 
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a/ A 
GGA 6.0279 


c/ A 
7.4631 


7 
120° 


V/A 3 
234.849 


EXP. 5.9037 


7.2636 


120° 


219.246 



Tabic 5: Lattice parameters for BCN compound calculated within GGA-PBE calculation 
using the norm-conserving pseudopotential. Experimental data (EXP.) for BCN ceramic 
are also presented. [25] 










UOOt 




750 800 850 

Wavenumber / cm' 1 



750 800 850 

Waveriumber / cm' 



900 



Figure 4: Deconvolution of the Raman spectra of the BCN ceramics in the spectral range 
700-900 cm -1 for several values of sintering temperature (a) and sintering time (b). Col- 
lected data (o) and calculated Lorentzian (-) curve by fitting. 
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Figure 5: Behavior of the ratio ^Ca,Nb in function of the sintering temperature and of the 
sintering time for BCN ceramics. 



sintering parameter changes. The most ordered sample is that sintered at 
1600 °C for 2h for which # Ca (0.03) and *m,(0.97). We observe the sintering 
time is not so efficient to order BCN, once that the calculated ratios to BCN 
ceramic sintered at 1300 °C for 16h were \Pc o (0.13) and ^1/^^(0.87). 

5. Conclusions 

Partially ordered BCN ceramics was obtained by polymeric precursor 
method and the ordering in these samples were investigated by Raman spec- 
troscopy under sintering temperature and sintering time changes. The evo- 
lution of the order was evaluated with basis on a partially ordered trigonal 
structure. An ab initio calculation permit us to assign the phonons and to 
monitor the order observing the changes in the behavior of the A\ g Raman- 
active mode near to 821 cm" 1 . The most ordered sample was obtained at a 
sintering temperature of 1600 °C at 2 h. 
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